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j^/ 1 Introduction 

^"J . The two archetypal ensembles of random matrices are Wigner real symmet- 

i-^h ' ric (Hermitian) random matrices and Wishart sample covariance real (corn- 

ed ■ plex) random matrices. In this paper we study the statistical properties of 

the largest eigenvalues of such matrices in the case when the second moments 
of matrix entries are infinite. In the first two subsections we consider Wigner 
ensemble of random matrices and its generalization - band random matrices. 

> 

1.1 Wiener Random Matrices 

in 

^■Lj ■ n x n matrix with i.i.d. entries up from the diagonal 

in ' 

. A — (ajk), o,jk — a kj, 1 < 3 ' < k < n, \Q>jk}j<k~ i-i.d. real random variables. 

(1) 

The diagonal entries {an}, 1 < i < n, are usually assumed to be i.i.d. random 
variables, independent from the off-diagonal entries. A Hermitian Wigner ran- 
Q dom matrix is defined in a similar way, namely as a square n x n Hermitian 

^ ' matrix with i.i.d. entries up from the diagonal 



A real symmetric Wigner random matrix is defined as a square symmetric 



■ A — (cijk), CLjk = a kj< 1 < j < k < n, {a,jk} j<k~ i-i.d. complex random variables. 

As in the real symmetric case, it is usually assumed that the diagonal entries 
{an}, 1 < i < n, are i.i.d. (real) random variables independent from the 
off-diagonal entries. 

Ensembles JQ and (J2J were introduced in mathematical physics by Eugene 
Wigner in the 1950s (03], ESI)- Wigner viewed these ensembles as a 

mathematical model to study the statistics of the excited energy levels of 
heavy nuclei. 
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The famous Wigner's semicircle law can be formulated as follows. Let the 
matrix entries in (JTJ or be centered random variables with the tail of 
distribution decaying sufficiently fast, so that all moments exist. Denote by 
Ai > A2 > . . . > A ra the eigenvalues of a random matrix rT 1 ! 2 A. Then the 
empirical distribution function of the eigenvalues converges, as n — > 00, to a 
non-random limit 

-#(A, < x, 1 < i < n) -> F(x) = f f{t)dt, (3) 



where the density of the semicircle law is given by fit) = ■^2\/2<j 2 — x 1 , for 
t G [— v2<7, v2cr]; and <r 2 is the second moment of matrix entries. 

This result was subsequently strengthened by many mathematicians (see 
e.g. pQ, [S], PS)- I n ^s general form (due to Pastur and Girko), the theo- 
rem holds if the matrix entries of A satisfy the Lindeberg-Feller condition: 
■h Ei<,<j<„ S\x\>T^i x2dF ij ( x ) ~* 0) where (x) is the distribution func- 
tion of af-i . 

From the analytical point of view, the simplest examples of Wigner ran- 
dom matrices are given by the so-called Gaussian Orthogonal and Unitary 
Ensembles (GOE and GUE for short). The GUE is defined as the ensemble of 
nxn Hermitian matrices with the Gaussian entries ~ N(0, 1 /2) , ^sajk ~ 
N(0, 1/2), 1 < j < k < n; a u ~ iV(0, 1), 1 < i < n (see ([37], chapter 6). 
The joint distribution of the matrix entries has the form 

P(dA) = const n exp ^—^Tr(A 2 )^j dA, (4) 

where dA = Yij<k d^Ra^dQajk Yii=i da>u is the Lebesgue measure on the space 
of n-dimensional Hermitian matrices. The joint distribution of the eigenvalues 
is given by its density 



Zn 1 II {xi-Xjfexp 

l<i<7*<n \ i / 



p n (xi, . . . ,x n ) = Z n [[ {xi -Xj) exp ( -- > ^ j . (5) 

i<j< 

The normalization constants in (jJJ and J^J) are known. What is more, one 
can calculate explicitely the /c-point correlation functions (see [27J> chapter 
6). This allows one to study the local distribution of the eigenvalues, both 
in the bulk of ths spectrum and at its edges in great detail. In particular, a 
celebrated result of Tracy and Widom (see 0U|) states that 

lirn^ Pr (x m ax < %Vn + ^TTe ) = F ^ = exp [~ f s ^ ~ s )^i x )^j > 



(6) 

where q(x) is the solution of the Painlcvc II differential equation 

q"(x) = xq{x) + 2q 3 {x) 
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with the asymptotics at infinity q(x) ~ Ai(x) at x — +00. 

The limiting /c-point correlation function at the edge of spectrum is given 
by the formula 

p k {xi,...,x k ) =det(K{x i ,x j )) 1 < itj < k , (7) 

where 

K(x, y) = K Airy (x, y) = M*W (y) - M> (x) M y) 

x-y 

is a so-called Airy kernel. We refer the reader to 0U| and J3| for the details. 
We recall that the fc-point correlation function is defined in such a way that 
for any disjoint subintervals of the real line I±, 1%, . . . , Ik, one has 

k 

E X\Wi) = l •••/ Pk(xi, ■ ■ ■ ,x k )dx 1 . . .dx k , 

where #(I) denotes the number of the eigenvalues in /. A probabilistic 
interpretation of the above formula is that p k {xi, . . . , x k )dxi . . . dx k is the 
probability to find an eigenvalue in each of the k infinitesimal intervals [x{ , Xi+ 
dxi], i=l,...,k. 

The Gaussian Orthogonal Ensemble (GOE) is defined as the ensemble 
ofnxii Wigner real symmetric random matrices with the Gaussian entries. 
More precisely, we assume that ay , 1 < i < j < n, are independent Gaussian 
iV(0, 1 + Sij) random variables (see e.g. chapter 7). The joint distribution 
of the matrix entries has the form 

P(dA) = c„ exp (~iTr(A 2 )^j dA, (9) 

where dA — Yii<j ^ a ij i- s the Lebesgue measure on the space of n-dimensional 
real symmetric matrices. The distribution (J5J) induces the joint distribution 
of the eigenvalues of the GOE matrix, given by its density 

p„(xi,...,x n ) =Z~ l Y[ \ x i ~ ^ilexp I -jX^) • ( 10 ) 

\<i<j<n \ i / 

The limiting distribution of the (normalized) largest eigenvalue of a GOE 
matrix was calculated by Tracy and Widom in (03). 

+00 

2/ 



(x) + (x — s)q (x)dx 



Jhn^Pr (\ max < 2^n + ^T/e) = = ex P J 1 

(11) 

The Tracy- Widom distribution pill was obtained by studying the asymptotic 
properties of the /c-point correlation functions at the edge of the spectrum. The 
fc— point correlation function in the GOE ensemble has the pfaffian form. In 
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the limit n — > oo the fc-point correlation function at the edge of the spectrum 
is given by the following formula 

S \ 1/2 
p k (xi,...,x k ) = [det{K(x l ,x j )) 1 ^ j ^ k J , (12) 

where K(x,y) is a 2 x 2 matrix kernel such that 

1 f y 

K n {x, y) = K 22 (y, x) = K Atry (x, y) + -Ai(x) / Ai(t)dt, (13) 

* J -oo 

1 d 

K 12 (x, y) = --Ai{x)Ai{y) - —K A „ y (x, y), (14) 

2 oy 

K 2 i(x,y) = / I / Ai{v)dv I Ai(x + u)du — e(x — y) + — / Ai(u)du 
Jo \Jx+u ) 2 J y 

+ - / Ai{u)du I Ai(v)dv, (15) 

2 Jx J -oo 

where e(z) = ^sign(z). 



1.2 Band Random Matrices 

A band random matrix is a generalization of a Wigner random matrix en- 
semble Q, (J2J. A real symmetric (aperiodic) band random matrix is defined 
as a square symmetric n x n matrix A = (ajk) such that et^ = unless 
\i - j\ < d„, and 

{ a jk-> j < k; \j — k\ < d n } — i.i.d. real random variables. (16) 

A Hcrmitian band random matrix is defined in a similar way, namely as 
a square n x n Hermitian matrix A = {a-jk), such that ay = unless 
\i — j\ < d n , and 

{ajfc: j < k; \j — k\ < d n } — i.i.d. complex random variables. (17) 

If d n = n — 1, we obtain the Wigner ensemble of random matrices. A 
matrix is called a periodic band matrix if \i — j\ is replaced above by \i — j\i = 
min(|i — j\,n — \i — j\). Band random matrices have been studied in the 
last fifteen years (see for example 0, 0, 123 > I n the periodic case, the 

limiting distribution of the eigenvalues of (i" 1 A is given by the semi-circle law, 
provided matrix entries have a finite second moment. In the aperiodic case, 
the limiting distribution of the eigenvalues is different from the semi-circle 
law, unless d n /n — ► (see e.g. [SOI)- One of the most interesting problems 
involving band random matrices is the localization/ derealization properties 
of the eigenvalues. It is conjectured in physical literature, that the eigenvalues 
of band random matrices are localized if d n — 0(n 1 ^ 2 ). As far as we know, 
there are no rigorous results yet in this direction. 
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1.3 Sample Covariance Random Matrices 

Sample covariance random matrices have been studied in mathematical statis- 
tics for the last seventy- five years. We refer to [SI and 53 f° r the ap- 
plications of spectral properties of Wishart random matrices in multivariate 
statistical analysis. 

Let A be a large m x n real rectangular random matrix with independent 
identically distributed entries. In applications, one is often interested in the 
statistical behavior of the singular values of A in the limit m — > oo, n — > oo. 
This is equivalent to studying the eigenvalues of a positive-definite matrix 
M = A 1 A in the limit of large dimensions. Without loss of generality, one 
can assume that m > n (since the spectrum of AA differs from the spectrum 
of A 1 A only by a zero eigenvalue of multiplicity m — n. 

The analogue of the Wigner semicircle law was proved by Marchenko and 
Pastur (|23)- Let m — ► oo, n — > oo in such a way that m/n — > 7 > 1. Assume 
E\a,ij\ 2+e < +00, where e > is an arbitrary positive number. Then the 
empirical distribution function of the eigenvalues of ^A t A converges to a 
non-random limit, known as the Marchcnko-Pastur distribution 

1 f x 

-#(Ai <x, i= l,...,n) ^G 7 (x) = / g 7 (t)dt, (18) 

n J -co 

where the spectral density g(t) is supported on the interval [a, b], a = a 2 (l — 

7 -V2)2 j b = ff 2 (1+7 -l/2 ) 2 j ff 2 = andff(t) = _L_^ v / (6 _ t)( i_ a); t g 

[a, 6]. 

The ~ AT(0, 1) 1 < i, j < n, is known in the literature as the 

Wishart (Laguerre) ensemble of real sample covariance matrices. The joint 
distribution of the eigenvalues of M is defined by its density. Similarly to the 
Gaussian ensembles of real symmetric and Hermitian matrices discussed in 
Subsection 1.1, many important statistical quantities in the Wishart ensemble 
can be calculated explicitely. For example, the joint probability density of the 
eigenvalues is given by the formula 

n 

p n (x u ...,x n ) = Z-i m J] Ixi-x^jlxr-^expi-Xi/i). (19) 

l<z<j<n i—1 

It was shown by Johnstone that the largest eigenvalue of a Wishart 

random matrix converges, after a proper rescaling, to the Tracy- Widom dis- 
tribution F\. Namely, let m — > 00, n — > 00, m/n — * 7 and u m n — 

( n l/2 +m l/2) 2) Umn= (n l/2 +m l/2 )(n -l/ 2 + m -l/2 ) l/3 > Then 

)^F 1 (s) (20) 

One can also show (see jSZ]), that the rescaled fc-point correlation function 
at the edge of the spectrum converge in the limit to (|12H . 

Finally, we want to remark, that there is a long-standing interest in nuclear 
physics in the spectral properties of the complex sample covariance matrices 
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A* A, where the entries of a reactangular matrix A are independent identically 
distributed complex random variables (see e.g. 013, [S], 132 ■ El)- We refer 
the reader to |37| and the references therein for additional information. 

1.4 Universality in Random Matrices 

The universality conjecture in Random Matrix Theory states, loosely speak- 
ing, that the local statistical properties of a few eigenvalues in the bulk or at 
the edge of the spectrum are independent of the distribution of individual ma- 
trix entries in the limit of large dimension. The only thing that should matter 
is, whether the matrix is real symmetric, Hermitian or self-dual quaternion 
Hermitian. 

For Wigner random matrices, the conjecture was rigorously proven at the 
edge of the spectrum, both for real symmetric and Hermitian case in |36|. 
provided that all moments of matrix entries exist and do not grow faster than 
the moments of a Gaussian distribution, and the odd moments vanish. In 
particular, it was shown that the largest eigenvalue, after proper rescaling, 
converges in distribution to the Tracy- Widom law. In the bulk of the spec- 
trum, the conjecture was proven by Johansson ( 20 ) for Wigner Hermitian 
matrices, provided the marginal distribution of a matrix entry has a Gaussian 
component. We refer to 1 1 1 j and references therein for the universality results 
in the unitary ensembles of random matrices. 

The situation for sample covariance random matrices is quite similar (see 
papers by Soshnikov [37j and Ben Arous and Pcchc 3 ). 

The natural question is how general such results are? What happens if ma- 
trix entries have only a finite number of moments? In this article we consider 
the extreme case when the entries of A do not have a finite second moment. 
In the next section, we discuss spectral properties of Wigner random matrices 
and, more generally, band random matrices when marginal distribution of ma- 
trix entries has heavy tails. As was shown in PHj , the statistics of the largest 
eigenvalues of such matrices are given by a Poisson inhomogeneous random 
point process. In Section 3 we discuss a similar result (although in a weaker 
form) for the largest eigenvalues of sample covariance random matrices with 
Cauchy entries. Section 4 is devoted to conclusions. 



2 Wigner and Band Random Matrices with Heavy Tails 
of Marginal Distributions 

In this section we consider ensembles of Wigner real symmetric and Hermitian 
matrices JQl and J2J > an d band real symmetric and Hermitian random matrices 
HI tip , (PJl with the additional condition on the tail of the marginal distribution 

G(x)=Pr(\a jk \>x) = ^, (21) 
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where < a < 2 and h(x) is a slowly varying function at infinity in a sense 
of Karamata ( E^J , EH)- I n other words, h(x) is a positive function for all 
x > 0, such that lim^oo = 1 for all t > 0. The condition (|21|l means 
that the distribution of |a,j| belongs to the domain of the attraction of a 
stable distribution with the index a (see e.g. JB], Theorem 2.6.1). 

Without loss of generality, we restrict our attention to the real symmet- 
ric case. The results in the Hermitian case are practically the same. Wigner 
random matrices JTJ, with the heavy tails 121|) . in the special case when 
limit lim^^oo h(x) > exists, were considered on a physical level of rigor 
by Cizeau and Bouchaud in They argued, that the typical eigenvalues of 
A are of the order of n x l a . Cizeau and Bouchaud also suggested a formula 
for the limiting spectral density of the empirical distribution function of the 
eigenvalues of n~ 1 / a A. Unlike the Wigner semicircle and Marchenko-Pastur 
laws, the conjectured limiting spectral density is supported on the whole real 
line. It was given as 

f{x) = L C Jf^\x), (22) 

where L^'^ is a density of a centered Levy stable distribution defined through 
its Fourier transform L{k) : 

L C J = tM dkL(ky k *, (23) 



2vr . 

lni(fc) = -C\k\ a (l + if3sgn(k)t&n(Tra/2)), (24) 
and functions C(x), j3[x) satisfy a system of integral equations 

C(x)= / \y\^- 2 L C a f^\x-^)dy, (25) 



-OO 



P(x)=l L C a ^\x--)dy. (26) 



We would like to draw the reader's attention to the fact that the density 
in (121'li is not a density of a Levy stable distribution, since C(x), 0(x) are 
functions of x. Cizeau and Bouchaud argued, that the density f(x) should 
decay as xi \ a at infinity, thus suggesting that the largest eigenvalues of A 
(in the case h(x) — const) should be of order ni , and not n«, which is the 
order of typical eigenvalues. 

Even though originally proven in |38| in the Wigner case, the theorem 
written below holds in the general case of band random real symmetric (or 
Hermitian) random matrices l|16|l . (|17|l . 

Let N n be the number of independent (i.e. i < j), non-zero matrix entries 
ay in A. In other words, let N n = =ff(l < i < j < n, \i — j\ < d n ) in 
the aperiodic band case, and iV„ = ^=(1 < i < j < n, \i — j\i < d n ) in 
the periodic band case. It is not difficult to see, that N n = i n the 
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Wigner case, N n = n(d n + 1) in the periodic band case, and N n = n x 
(rin + 1)- d " (d 2 " +1) in the aperiodic band case. Let us define a normalization 
constant b n in such a way that 

lim N n G{b n x) = ^-, (27) 

n — >oo X 

for all positive x > 0, where the tail distribution G has been defined in (|21|1 . 
Normalization b n naturally appears (see [231 an d Remark 1 below), when 
one studies the extremal values of a sequence of N n independent identically 
distributed random variables i|21|) . In particular, one can choose 

b n = inf{t : G(i - 0) > > G(t + 0). (28) 

It follows from (23 and that iV^" 5 < 6„ < A^ +<5 for arbitrary small 
positive 5, and Nnl ^ b ' l " > — > i as n — > oo. 

Theorem 1 claims that the largest eigenvalues of A have Poisson statistics 
in the limit n — > oo. 

Theorem 1. Let A be a band real symmetric if_?6)) or Hermitian j_?7| ) ran- 
dom matrix with a heavy tail of the distribution of matrix entries Then 
the random point configuration composed of the positive eigenvalues of b^A 
converges in distribution on the cylinder sets to the inhomogeneous Poisson 
random point process on (0, +oo) with the intensity p(x) — i" a . 

In other words, let < x\ < y\ < x-i < yi < . . . xj- < yk < +oo, and 
Ij = i x 3' Ui)i J = 1 , . . . , be disjoint intervals on the positive half-line. Then 
the counting random variables = #(1 < i < n : A; G Ij), j = 1, . . . , k, 

are independent in the limit n — > oo, and have a joint Poisson distribution 
with the parameters fij — jj p(x)dx, i.e. 

fe Sj 

lim Pr(#(/ i )= ail j = l,...,A) = TT^ r e-w. (29) 

i=i J 

For the additional information on Poisson random point processes we refer 
the reader to [TU|. 

Corollary 1. Let be the k-th largest eigenvalue ofb~ l A, then 

fe— l —i a 

lim Pr(A fe < x) = exp(-2T Q ) V ^— -. (30) 

Z=0 

In particular, lim„_ > . 0O Pr(Ai < x) — exp(— x~ a ). 
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Remark 1 The equivalent formulation of the theorem is the following. Let 
k be a finite positive integer. Then the joint distribution of the first k largest 
eigenvalues ofb^A is asymptotically (in the limit n — » oo) the same as the 
joint distribution of the first k order statistics of {b^ |ffly|, 1 < i < J < ra}. 
It is a classical result, that extremal values of the sequence of independent 
identically distributed random variables with heavy tails distributions $21\) have 
Poisson statistics (see e.g. [23], Theorem 2.3.1). 

Theorem 1 was proven in jSH] in the Wigner (i.e. full matrix) case ©■ 
The proof in the general (band matrix) case is essentially the same. However, 
it should be noted, that the original proof of Theorem 1 in contained a 
little mistake, which could be easily corrected. The corrections are due in two 
places. 

First of all, the correct formulation of the part c) of Lemma 4 from 38 (p. 
87) should state, that for any positive constant 5 > 0, with probability going 
to 1 there is no no row 1 < i < n, that contains at least two entries greater 

in absolute value than b„ +S ■ In other words, the exponent \ + 5 in bn +S 
in part c) of Lemma 4 must be replaced by | + 6. After this correction, the 
statement is true. Indeed, the probability that there is a row with at least two 

entries greater than b r \ +S can be estimated from above by n 3 (^G(bn +S )j . 

It follows from (|21[1 . (|27[1 and i|28|) . that this probability goes to zero. 
Also, the formula (28) in Lemma 5 (p. 88) should read 

Pr{3i, 1 < i < n : max |<2j,-| > ( |ay|)— max \a i3 \ > b^ +T } — > 

l<j<n ^ — ' l<j<n 
l<j<n 

(31) 

as n — > oo. In other words, the exponent i + j in b% + 4 must be replaced by 
I + f ■ The key step of the proof of Lemma 5 was to show, that for any fixed 
row i and arbitrary small positive 5, the probability Pr(J^ i +s |oy| > 

3':|o«l<6n 

bn +2 ) can be estimated from above by exp(— n e ), where e = e(S, a) > 0. We 
then concluded, that with probability going to 1, there is no row i such that 

i+« \ a ij\ — b?i ■ To establish (|3*Hl . it is enough to prove that for 

any fixed row i 

Pr( V \(Hj\ >bl +2& ) <exp(-n £ ), (32) 

i:by 2+5 <|a 13 |<6| + 5 

for sufficiently small positive e. The proof is very similar to the argument 
presented in Lemma 5, and is left to the reader. 

3 Real Sample Covariance Matrices with Cauchy Entries 

Let A be a rectangular mxn matrix with independent identically distributed 
entries with the marginal probability distribution of matrix entries satisfying 
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J2U- Based on the results in the last section, one can expect that the largest 
eigenvalues have Poisson statistics as well. At this point, we have been able 
to prove it only in a weak form, and only when matrix entries have Cauchy 
distribution. 

We recall, that the probability density of the Cauchy distribution is given 
by the formula f(x) — -^j-pjry- Cauchy distribution is very important in 
probability theory (see e.g. 53)- In particular, Cauchy distribition is a (1,1,0) 
stable distribution, i.e. the scale parameter is 1, the index of the distribution 
a = 1 and the symmetry parameter is zero (see |18j . |23jV 

The following theorem was proven by Fyodorov and Soshnikov in |39) 

Theorem 2. Let A be a random rectangular to x n matrix (m > n ) with 
i.i.d. Cauchy entries and let z be a complex number with a positive real part. 
Then, as n — > oo we have 

^ 1/2 f o \ 00 

lim E (det(l H ^—zA t Ay\ = exp ( ^ ) = ETT(1 + zx,)- 1/2 , 

v ' i=l 

(33) 

where we consider the branch of ^pz on D — {z : Sftz > 0} such that vT = 1, 
E denotes the mathematical expectation with respect to the random matrix 
ensemble defined above, E denotes the mathematical expectation with respect 
to the inhomogeneous Poisson random point process on the positive half-axis 
with the intensity 3/2 , an d the convergence is uniform inside D (i.e. it is 
unform on compact subsets of D). For a real positive z — t 2 , t € R , one can 
estimate the rate of convergence, namely 

lim E (&ct(l H i ^—A t A)] = cxp ( -~\t\(l + ofa" 1 ^)) V (34) 

where e is an arbitrary small positive number and the convergence is uniform 
on the compact subsets of [0, +00). 

The result of Theorem 2 allows a generalization to the case of a sparse 
random matrix with Cauchy entries. Let, as before, {ajfe}, 1 < j < m, 1 < 
k < n, be i.i.d. Cauchy random variables, and B = (bjk) be a m x n non- 
random rectangular — 1 matrix such that the number of non-zero entries in 
each column is fixed and equals to d n . Let d n grow polynomially, i.e. b n > n a , 
for some < a < 1. Also assume that lnm grows much slower than than 
any power of n. 

Define a m x n rectangular matrix A with the entries Ajk — bjkOjk, 1 < 
j < to, 1 < fc < n. Let Ai > A2 . . . > A„ denote the eigenvalues of A 1 A. The 
appropriate rescaling for the largest eigenvalues in this case is Ai = — ^ , i — 
l,...,n. 

Theorem 3. Let A be a sparse random rectangular to x n matrix (m > n) 
defined as above, and let z be a complex number with a positive real part. 
Then, as n — > 00 we have 
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\ -1/2 n 
lira E ( det(l H ^- Fr A t A) ) = lira J^TTfl + zAj)~ 1/2 (35) 

n. — kyi \ mAnA 1 n — >oo J- J- 



m 2 (R, ' I n— >oo 

" ' i=l 



2 x 



= exp I —^l J = E JJ(1 + z^)" 1 / 2 , (36) 
^ 77 ' i=l 

where, as in Theorem 1.1, we consider the branch of -^/z on D — {z : -Rz > 0} 
smc/i i/ia£ \/T = 1. E denotes the mathematical expectation with respect to the 
inhomogeneous Poisson random point process on the positive half-axis with 
the intensity —^j? , and the convergence is uniform inside D (i.e. it is unform 
on the compact subsets of D). For a real positive z = t 2 , t& R , one can get 
an estimate on the rate of convergence, namely 

e ( det(i + ^ n AtA) ) 1/2 = exp H <(i + °^ i/2+£ ))) • <w 

where e is an arbitrary small positive number and the convergence is uniform 
on the compact subsets of [0, +00). 

The proof relies oa the following property of the Gaussiaa iategrals: 

(det(B)) _1/2 = Q J xexp(-xBx t )d N , (38) 

where B is an TV-dimensional matrix with a positive definite Hermitian part 
(i.e. all eigenvalues of B + B* are positive), x — (x\, . . . ,xn) & R N , and 

xBx — y ] jj bijXiXj . 

Let B = B(t) = (fyj^t ^jd^J ' Then, one can write 

(det(l + t 2 ^))- 1/2 = (det^*^ 4 )) ^^(det^))- 172 , (39) 

and apply l|38|l to the r.h.s. of (|39|l . Assumiag that the eatries of A are 
iadepeadeat, oae caa sigaificaatly simplify the expression, using the fact 
that entries of A appear liaearly ia B(t) (see Proposition 1 of [SI). Ia the 
Cauchy case, oae caa simplify the calculatioas evea further, aad prove that 

lioin^oo E ( det(l + -^§ 1 rA t A)J exists and equals exp (-~y/z) . 

On the other side, for Poisson random point processes the mathematical 
expectations of the type Ef]™^ + f(%i)) can be calculated explicitely 

00 00 k 

Ena + /^)) = 1 +E E E II = ( 4 °) 

i=l k—1 l<ii<i2<-..<ik 3 — 1 

00 2 /■ 

= EtT/ J] f(xj)pk(xi,...,Xk)dxi---dx k (41) 

k=0 k] - ./(o.+oo)* fJi 
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= E h ( I f( x )p( x ) dx ) = exp f / f(x) P (x)dx) . (42) 

fc=0 \- / (0,+<») / \->(0,+°o) J 

In the equations above, pk denotes the /c-point correlation function, and p 
denotes the one-point correlation function (also known as intensity). It is a 
characteristic property of a Poisson random point process that the k-point 
correlation function factorizes as a product of one-point correlation functions, 
i.e. pk(xi, ■ ■ ■ ,Xk) = rii=i P( x i)- I 11 the context of Theorems 2 and 3, test 
function / has the form f(x) = (1 + zx)^ 1 ^ 2 — 1. When the intensity p equals 
— i-ra, one obtains 

{ f(x)p(x)dx= { {{l + zx)- l/2 -\)—Lj-dx = --y/z, 

J(Q,+oc) J(0,+oo) 6 I* 7T 

which finishes the proof. 

The fact, that the intensity p[x) = — of the Poisson random point 
process diverges at zero and is summable at +00, means, that the the vast 
majority of the eigenvalues of the normalized matrix converge to zero in the 
limit. 

Remark 2 It should be pointed out, that the results of Theorem 2 and 3 do 
not imply that the statistics of the largest eigenvalues of a normalized sample 
covariance matrix with Gaussian entries are Poisson in the limit of n — > 00. 
Indeed, to prove the Poisson statistics in the limit one has to show that 

n +00 

n lim ^n( 1 + /(^)) =E II( 1 + ^)) ( 43 ) 

i=l i=l 

for a sufficiently large class of the test functions f , e.g. for step functions with 
compact support. As we already pointed out, the results of Theorems 2 and 3 
claim that 14 3[ ) is valid for f(x) — (1 + zx)^ 1 / 2 — 1 for all z such that 5Rz > 0. 



4 Conclusion 

It is known in the theory of random Schrodinger operators, that the statistics 
of the eigenvalues is Poisson in the localization regime (see e.g. j2Hl> EH1)- It 
seems, that the same mechanism is responsible for the Poisson statistics for 
the largest eigenvalues in the random matrix models described above. The 
interesting next problem is to find a phase transition between the Tracy- 
Widom regime (when all moments of matrix entries exist) and the Poisson 
regime (when second moment does not exist). 

It is also worth to point out, that there is a vast literature on the Poisson 
statistics of the energy levels of quantum sysytems in the case of the regular 
underlying dynamics (see e.g. 0], [H], 
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